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The possibility of having an extension of the cosmological standard model with a Hubble expan- 
sion rate H constrained to a finite interval is considered. Two periods of accelerated expansion arise 
naturally when the Hubble expansion rate approaches to the two limiting values. The new descrip- 
tion of the history of the universe is confronted with cosmological data and with several theoretical 
ideas going beyond the standard cosmological model. 



INTRODUCTION 

According to General Relativity, if the universe is filled 
with the particles of the Standard Model of particle 
physics, gravity should lead to a deceleration of the ex- 
pansion of the universe. However, in 1998 two indepen- 
dent evidences of present accelerated expansion were pre- 
sented [TJ [5] and later confirmed by different observations 
[3J|31[S]. On the other hand, measurements of large scale 
structure [6] and CMB anisotropy [3] also indicate that 
the universe evolved through a period of early accelerated 
expansion (inflation). 

There is no compelling explanation for any of these 
cosmic accelerations, but many intriguing ideas are being 
explored. In the case of inflation, the origin of the accel- 
erated expansion can be either a modification of gravity 
at small scales [7] or a coupling of the expansion of the 
universe to the progress of phase transitions [H] . In the 
case of the present accelerated expansion these ideas can 
be classified into three main groups: new exotic sources 
of the gravitational field with large negative pressure [5] 
(Dark Energy), modifications of gravity at large scales 
[TU] and rejection of the spatial homogeneity as a good 
approximation in the description of the present universe 

nu. 

Different models (none of them compelling) of the 
source responsible for each of the two periods of acceler- 
ated expansion have been considered. Einstein equations 
admit a cosmological constant A, which can be realized as 
the stress-energy tensor of empty space. This A, together 
with Cold Dark Matter, Standard Model particles and 
General Relativity, form the current cosmological model 
(ACDM). However, quantum field theory predicts a value 
of A which is 120 orders of magnitude higher than ob- 
served. Supersymmetry can lower this value 60 orders of 
magnitude, which is still ridiculous [T^]. In order to solve 
this paradox, dynamical Dark Energy models have been 
proposed. 

This has also lead to explore the possibility that cosmic 
acceleration arises from new gravitational physics. Also 
here several alternative modifications of the Einstein- 
Hilbcrt action at large and small curvatures [13j [HI 
HSJ [TU [17l HBJ [19], or even higher dimensional models 
[2TH [2Tj , producing an accelerated expansion have been 



identified. All these analysis include an ad hoc restriction 
to actions involving simple functions of the scalar curva- 
ture and/or the Gauss-Bonnet tensor. This discussion is 
sufficient to establish the point that cosmic acceleration 
can be made compatible with a standard source for the 
gravitational field, but it is convenient to consider a more 
general framework in order to make a systematic analy- 
sis of the cosmological effects of a modification of general 
relativity. 

In this paper we parametrize the evolution of the uni- 
verse (considered isotropic and homogeneous) with the 
Hubble parameter H. One finds that it is possible to 
restrict the domain of H to a bounded interval. This 
restriction naturally produces an accelerated expansion 
when the Hubble expansion rate approaches any of the 
two edges of the interval. Therefore, we find a new way 
to incorporate two periods of cosmic acceleration pro- 
duced by a modification of general relativity. But the 
dependence on the two limiting values of H can be cho- 
sen independently. One could even consider a Hubble ex- 
pansion rate constrained to a semi-infinite interval with 
a unique period of accelerated expansion. In this sense 
the aim to have a unified explanation of both periods 
of accelerated expansion is only partially achieved. This 
simple phenomenological approach to the problem of ac- 
celerated expansion in cosmology proves to be equivalent 
at the homogeneous level to other descriptions based on 
modifications of the Einstcin-Hilbert action or the intro- 
duction of exotic components in the matter Lagrangian. 

In the next section we will review how a general modi- 
fication of the gravitational action leads to a generalized 
first Friedman equation at the homogeneous level. In 
the third section we will present a specific model (ACM) 
based on the simplest way to implement a bounded in- 
terval of H. In the fourth section we will contrast the 
predictions of ACM for the present acceleration with as- 
trophysical observations. In the fifth section we will show 
that it is always possible to find modified gravitational 
actions which lead to a given generalized first Friedman 
equation, and present some simple examples. In the sixth 
section we will show that it is also always possible to find 
Dark Energy models which are equivalent to a given gen- 
eralized first Friedman equation, and present some simple 
examples. The last section is devoted to summary and 
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conclusions. 



ACTION OF THE COSMOLOGICAL STANDARD 
MODEL EXTENSION 

The spatial homogeneity and isotropy allow to reduce 
the gravitational system to a mechanical system with two 
variables a(i), N(t) which parametrize the Robertson- 
Walker geometry 



ds 2 



N 2 {t)dt 2 - a 2 (t)A ij dx i dx j 



1 - fcx 2 



(1) 



(2) 



Invariance under parameterizations of the time variable 
imply that the invariant time differential N(t)dt must 
be used. Also a rescaling of the spatial variables x % — > 
\x l together with a(t) — > A _1 a(i) and k — ► A~ 2 fc is a 
symmetry that must be kept in the Lagrangian. The 
action of the reduced homogeneous gravitational system 
can be then written as 



I g = I dtNL 



with H = 4f^y- 

aN at 



k 



1 dH 1 d 

N~dT'Ndt 



1 dH 
N~dt 



(3) 



If we keep the standard definition of the gravitational 
coupling and the density (p) and pressure (p) of a cosmo- 
logical homogeneous and isotropic fluid as a source of the 
gravitational field, we have the equations of the reduced 
system 
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P = 



iTrG N ) p = - 



1 SIg 

a 3 SN(t) 



1 SI g 
No? 5a(t) 



(4) 



(5) 



It is possible to choose a new time coordinate t such that 



dt 

If 



(6) 



This is equivalent to set N — 1 in the action (|3| of the 
gravitational system and the evolution equations of the 
cosmological model reduce to a set of equations for the 
scale factor a(t). If we introduce the notation 



dt 



H 



(SiL) 



U) _ 



[a 3 d t L] , 



(7) 



(8) 



with d L L denoting the partial derivatives of 
the Lagrangian as a function of the variables 
(k/a 2 ,H,HM,H( 2 \...) we have 



/ 87rGjv 



oo i—2 
i=2 j=0 



(-iy 



-H^ (StL) 



(i-J-2) 



(9) 



(8ttG n )p 



i=2 



(10) 

In the homogeneous and isotropic approximation, the 
vanishing of the covariant divergence of the energy- 
momentum tensor leads to the continuity equation 

d , q, d 



dt 



[pa 



~ P dt a 



3 



In the radiation dominated era one has 
P = 3p 



1a 



(11) 



(12) 



where a is a constant parameterizing the general solution 
of the continuity equation. In a period dominated by 
matter one has a pressureless fluid and then 



p = 0. 



(13) 



with constant rj. When these expressions for the energy 
density and pressure are plugged in (9p0l, one ends up 
with two compatible differential equations for the scale 
factor a(t) which describe the evolution of the universe. 
From now on we will use the more common notation 

dH — u 
St = 

THE ASYMPTOTIC COSMOLOGICAL MODEL 

Let us assume that the Lagrangian L of the gravita- 
tional system is such that the evolution equations (9fl0 1 
admit a solution such that a(t) > (absence of singu- 
larities), a > (perpetual expansion) and H < 0. In 
that case, one has a different value of the scale factor a 
and the Hubble rate H at each time and then one has a 
one to one correspondence between the scale factor and 
the Hubble rate. Since the continuity equation (together 
with the equation of state) gives a relation between the 
scale factor and the density, one can describe a solution 
of the evolution equations of the generalized cosmological 
model through a relation between the energy density and 
the Hubble rate, i.e. through a generalized first Friedman 
equation 



8ttG 



9(H) 



(14) 



with g a smooth function which parametrizes the dif- 
ferent algebraic relations corresponding to different so- 
lutions of different cosmological models. Each choice for 
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the function g(H) defines a phenomenological description 
of a cosmological model. Then one can take it as a start- 
ing point trying to translate any observation into a par- 
tial information on the function g(H) which parametrizes 
the cosmological model. 

Eq. (14 1 is all one needs in order to reconstruct the 



evolution of the universe at the homogeneous level. The 
generalized second Friedman equation is obtained by us- 



ing the continuity equation (111 and the expression for p 
as a function of H . One has 



inG N )p = 3g(H) + 9 -^P-H 



H 



In the matter dominated era one has 



H 
IP 



= -3 



9(H) 
Hg'(H) 



(15) 



(16) 



and then the assumed properties of the solution for the 
evolution equations require the consistency conditions 



g(H) >0, 



g (H) > . 



In the period dominated by radiation one has 



H 
H 2 



, 9(H) 
Hg'(H) 



(17) 



(18) 



instead of (16) and the same consistency conditions (17) 



for the function g(H) which defines the generalized first 
Friedman equation. 

We introduce now a phenomenological cosmological 
model defined by the condition that the Hubble rate has 
an upper bound H + and a lower bound This can 

be implemented through a function g(H) going to infin- 
ity when H approaches H + and going to zero when H 
approaches ii_ . We will also assume that there is an in- 
terval of H in which the behavior of the energy density 
with the Hubble parameter is, to a good approximation, 
scale- free i.e. g(H) ex H 2 . The source of the gravitational 
field will be a homogeneous and isotropic fluid composed 
of relativistic and non-relativistic particles. The Cosmo- 
logical Standard Model without curvature is recovered in 
the limit H- /H — -» and H + /H — » oo which is a good 
approximation for the period of decelerated expansion. 

Notice that this interpretation is independent of the 
underlying theory of gravitation. The Hubble parameter 
can be used to parametrize the history of universe as long 
as H ^ 0V£. The total density can be thus expressed as 
a function of H . If H is bounded, then p(H) will have 
a pole at H = H + and a zero at H = H-. Far from 
these scales, the behavior of p(H) can be assumed to be 
approximately scale-free. 

Under these conditions, we can parametrize the depen- 
dence of the cosmological model on the lower bound H_ 

by 



l(H) = H 2 h. 



H 2 _ 
IP 



(19) 



and similarly for the dependence on the upper bound H 

H 2 



9(H) = 



where the two functions h± satisfy the conditions 



(20) 



lim h±(x) = [3 

x^Q 



±1 



lim h±(x) = 0. (21) 

X— »1 



(3 is a constant allowed in principle by dimensional ar- 
guments. If (3 1 then it can be moved to the lhs of 
the Friedman equation, turning Gn ^ (3Gm, and can 
be interpreted as the ratio between an effective cosmo- 
logical value of the gravitational coupling and the value 
measured with local tests. But (3 ^ 1 would be in con- 
flict with Nucleosynthesis, through the relic abundances 
of 4 He and other heavy elements (for 3 neutrino species) 
[22] , so we set 0=1. Therefore 

limh±(x) = l, ]imh ± (x) = 0. (22) 

x — >0 x — >1 

The consistency conditions (|17|) result in 



h±(x) > , h ± (x) > xh ± (x) (23) 

for the two functions h± defined in the interval < x < 1 . 
Thus, we can divide the cosmic evolution history into 
three periods. In the earliest, relativistic particles domi- 
nate the energy density of the universe and the general- 
ized first Friedman equation shows a dependence on the 
upper bound H+. There is also a transition period in 
which the effect of the bounds can be neglected and the 
rhs of the first Friedman equation is scale- free; this pe- 
riod includes the transition from a radiation dominated 
universe to a matter dominated universe. In the third 
present period, non-relativistic particles dominate the en- 
ergy density of the universe but the dependence on the 
lower bound H— must be accounted for in the generalized 
first Friedman equation. 

From the definition of the Hubble parameter, one has 



H 
IP 



aa 



(24) 



Then, in order to see if there is an accelerated or decel- 
erated expansion, one has to determine whether H/H 2 
is greater or smaller than — 1. 

In the period dominated by radiation one has 



H 
IP 



1 



1 - 



xh , (x) 
h + (x) 



(25) 



with x = H 2 /H+, where we have used (20) assuming 



that only the dependence on the upper bound (H + ) of 
the Hubble parameter is relevant. A very simple choice 
for this dependence is given by 



9(H) 



H 2 



Ml 

Hi 



(26) 
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with o+ a (positive) exponent which parametrizes the 
departure from the standard cosmological model when 
the Hubble rate approaches its upper bound. With this 
simple choice one has a transition from an accelerated 
expansion for H 2 > + a+) into a decelerated ex- 

pansion when H 2 < H+/(l + aq_), which includes the 
domain of validity of the standard cosmological model 
(#«#+). 

In the period dominated by matter (which corresponds 
to lower values of the Hubble rate) we assume that only 
the dependence on the lower bound of the Hubble 

rate is relevant. Then one has 



H 
H 2 



xh_ (x) 
h-(x) 



with x = H 2 /H 2 . We can also consider a dependence on 
H— parametrized simply by an exponent a_ 



9(H) 



H 2 



H 2 _ 
H 2 



(28) 



With this choice one has a transition from a decelerated 
expansion for H 2 > i/ 2 (1 + a_/2) leaving the domain 
of validity of the standard cosmological model and enter- 
ing into an accelerated expansion when the Hubble rate 
approaches its lower bound for H 2 < H^_(l + a_/2). 

The possibility of describing p as a function of H is in- 
dependent of the existence of spatial curvature k. How- 
ever, in the kinematics of observables in the expanding 
universe we do need to specify the value of A;. In the 
rest of the paper we will assume that the universe is flat 
(k = 0), although the same analysis could be done for 
arbitrary k. 

The properties of the expansion obtained in this simple 
example (a period of decelerated expansion separating 
two periods of accelerated expansion) are general to the 
class of phenomenological models with a generalized first 



Friedman equation (14 1 with g(H) satisfying the consis- 



tency conditions ( 17 1 and a Hubble rate constrained to a 
finite interval. The specific part of the example defined 



by ( 26|28 l is the simple dependence on the Hubble rate 
bounds and the values of the Hubble rate at the transi- 
tions between the three periods of expansion. From now 
on we will name this description the Asymptotic Cos- 
mological Model (ACM). With respect to the late accel- 
erated expansion, ACM can be seen as a generalization 
of ACDM, which can be recovered by setting a_ = 1. 
It also includes an early period of exponential expansion 
which can be seen as a phenomenological description of 
the evolution of the universe at inflation in the homoge- 
neous approximation. 



Horizon problem 

One can see that the horizon of a radiation-dominated 
universe can be made arbitrarily large as a consequence of 
an upper bound on the Hubble parameter and in this way 
one can understand the observed isotropy of the cosmic 
microwave background at large angular scales. 

Let us consider the effect of the modification of the 
cosmological model on the calculation of the distance 
dh{tf,ti) of a source of a light signal emitted at time 
ti and observed at time tf 



d h {tf,ti) = ca(tf) 



dt 
W) 



(27) We have 



dt 



da 



~ 3 



-1/4 



9 (H)dH 
AHg(H) 3 / 4 ' 



(29) 



(30) 



where in the first step we have used the definition of the 
Hubble expansion rate H and in the second step we have 
used the relation between the scale factor a and H as 
given by Oil}. We are considering both times U and 



tt in the radiation dominated period. 
The distance dh is then given by 



dh(Hf,Hi) 



1 



Hi 



dHg (H) 



(31) 



If Hi is very close to H + (i.e. if we choose the time ti 
when the light signal is emitted well inside the period of 
accelerated expansion) then the integral is dominated by 
the region around Hi which is very close to H+. Then 
one can approximate in the integrand 



9(H) 



H 2 + 



1 



Ml 

H 2 



(32) 



On the other hand if the observation is made at a time 
tf within the domain of validity of the cosmological stan- 
dard model {H 2 < Hj < H%) then the factor g(Hf) 1 ^ 2 
in front of the integral can be approximated by Hf and 
then one has 



dh(Hf, Hi) 



=5^(1/2,-0/4), 



(33) 



where B z (m,n) is the incomplete Beta function, and it 
can be made arbitrarily large by choosing Hi sufficiently 
close to H + . In this way we see that a cosmological model 
with a finite interval of variation for H solves the horizon 
problem. 



CONSTRAINTS OF ACM BY OBSERVATIONS 

In this section we will carry out a more technical anal- 
ysis about how the astrophysical observations constrain 
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the parameter space of ACM in the matter dominance 
period. This analysis is based on the use of (assumed) 
standard candles, basically Type la Supernovae [23] and 
CMB [25 . These observations constrain the parame- 
ter space to confidence regions in which the combina- 
tion f2 m = (1 — H 2 ^ / H$) a ~ is constrained to be around 
one quarter. The consideration of both Type la SNe 
and CMB together favor a_ > 1.5. The results of this 
analysis can be seen in figures (FIG. 1-3). A reader not 
interested in technical details might well skip this section. 

We center our discussion of the experimental tests of 
the Asymptotic Cosmological Model in the late acceler- 
ated expansion produced when H approaches its lower 
bound H-. The vast amount of supernovae data col- 
lected by [23] and the data from the SDSS Baryon 
Acoustic Oscillation [23], the mismatch between total 
energy density and total matter energy density seen at 
CMB anisotropics [25] and the measurements of present 
local mass density by 2dF and SDSS [5] [27] compared 
with the measurements of Hp from the HST Cepheids 
[2"8] show that the universe undergoes a surprising accel- 
erated expansion at the present time. 

We will firstly confront the model with the Supernovae 
la data from Riess et al. and SNLS collaboration. The 
usefulness of the Supernovae data as a test of Dark En- 
ergy models relies on the assumption that Type la SNe 
behave as standard candles, i.e., they have a well de- 
fined environment-independent luminosity C and spec- 
trum. Therefore, we can use measured bolomctric flux 
J~ = 4~ra an d frequency to determine luminosity dis- 
tance cLl and rcdshift z. The luminosity distance is given 
now by 



ACM with dataset gold 



d L (z) = c(l 

c(l+z) fH(z) 



g'{H)dH 



dz 



(34) 



The computed value must be compared with the one 
obtained experimentally from the measured extinction- 
corrected distance moduli (hq — hlog\p{ ^ ) + 25) for 
each SN. The SNe data have been compiled in refer- 
ences [23]) and we have limited the lowest redshift at 
cz < 7000fcm/s in order to avoid a possible "Hubble 
Bubble" |2SJ |3D] . Therefore our sample consists on 182 
SNe. We will determine the likelihood of the parameters 
from a x 2 statistic, 



X 2 (H ,H. 



) - Hoa) 2 



(35) 

where a v is the dispersion in supernova redshift due to 
peculiar velocities (we adopt || v p ||= 400fcm/s in units 
of distance moduli) , a^ i is the uncertainty in the in- 
dividual measured distance moduli /io,i, and fi Pt i is the 
value of fig at computed with a certain value of the set 
of parameters Hq, H-,ct-. This % 2 has been marginal- 
ized over the nuisance parameter Hq using the adaptive 




ACM with dataset gold + Wro = 0.28 +/- 0.03 + RCtIB = 1.70 +/- 0.03 




FIG. 1: Confidence regions in parameter space of the Asymp- 
totic Cosmological Model (ACM) with (bottom) and with- 
out (top) priors from CMB and estimations of present matter 
energy density at la (green), 2a (red) and 3a (blue). The 
ACDM is inside the la region if we consider no priors, but 
is moved outside the 2a region when we confront the SNe 
data with CMB and estimations from present matter energy 
density. In this case a value a_ > 1.5 is favored. 



method in reference [21] ■ The resulting likelihood distri- 
bution function e~ x ' has been explored using Monte 
Carlo Markov Chains. We get a best fit of ACM at 



a_ = 0.36 and ^ = 0.95, for which X = 157.7. In 
contrast, fixing a_ = 1, we get il\ = ^tt — 0.66 and 

X 2 = 159.1 for the best fit ACDM. The confidence re- 
gions are shown in Fig. 1 (top). 

We can add new constraints for the model coming from 
measurements of the present local matter energy den- 
sity from the combination of 2dF and SDSS with HST 
Cepheids, rendering fi m = 0.28 ± 0.03; and the distance 
to the last scattering surface from WMAP, which leads 
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Best fit ACM vs. LCDM uith priors 




FIG. 2: Distance Moduli vs. Redshift comparison between 
the best fits of ACM (H = 66 Km/ sMPc, red) and ACDM 
(Ho — 65 Km/sAlPc, blue) to the SNe data with priors. 
ACM fits clearly better the medium redshift SNe. 



ACM with dataset gold + RCMB = 1.70 +/- 0.03 




FIG. 3: Q, m vs. a confidence regions with the CMB prior only. 
The confidence regions show tt m = 0.25 ± 0.03 and a > 1.15 
at 2cr level. 



to r C MB = Vn m fpj = 1.70 ± 0.03 [23]. Includ- 
ing these priors we get a best fit of ACM at a_ = 3.54 
and -p=^ = 0.56 (our simulation explored the region with 
a_ < 3.6), for which \ 2 — 164.1. In contrast, we get 
tt A = 0.73 and X 2 = 169.5 for the best fit ACDM, which 
is outside the 2a confidence region shown in Fig. 1 (bot- 
tom). The fits of the best fit ACDM and ACM taking 
into account the priors to the SNe data are compared in 
Fig. 2. The information which can be extracted from the 
data is limited. This can be seen in Fig. 3, in which it is 
explicit that the data constrain mainly the value of the 

H 2 

present matter energy density, fi m = ^ = (1 — -fpt) a ~ ■ 



GENERALIZED FIRST FRIEDMAN EQUATION 
AND f(R) GRAVITY 

The extension of the cosmological model considered 
here could be compared with recent works on a mod- 
ification of gravity at large or very short distances 

naming nana Hang, it has been shown that, b y 

considering a correction to the Einstein Hilbert action in- 
cluding positive and negative powers of the scalar curva- 
ture, it is possible to reproduce an accelerated expansion 
at large and small values of the curvature in the cosmo- 
logical model. Some difficulties to make these modifica- 
tions of the gravitational action compatible with the solar 
system tests of general relativity have lead to consider a 
more general gravitational action, including the possible 
scalars that one can construct with the Riemann cur- 
vature tensor [32], although the Gauss-Bonnet scalar is 
the only combination which is free from ghosts and other 
pathologies. In fact, there is no clear reason to restrict 
the extension of general relativity in this way. Once one 
goes beyond the derivative expansion, one should con- 
sider scalars that can be constructed with more than two 
derivatives of the metric and then one does not have a 
good justification to restrict in this way the modification 
of the gravitational theory. It does not seem difficult to 
find an appropriate function of the scalar curvature or the 
Gauss-Bonnet scalar, which leads to a cosmology with a 
bounded Hubble expansion rate. 

One may ask if a set of metric f(R) theories which 
include ACM as an homogeneous and isotropic solution 
exists. The answer is that a bi-parametric family of f(R) 
actions which lead to an ACM solution exists. In the 
following section we will derive them and we will discuss 
some examples. The derivation follows the same steps 
of modified f(R)-gravity reconstruction from any FRW 
cosmology [3"3"] . 

We start from the action 

S = f d i x^[f(R) + l^GL m ] , (36) 

167tGat J 

which leads to the "generalized first Friedman equation" 



-18(HH + 4H 2 H)f"(R) 
-3(H + Hi)f>(R)-y(R) = (l 

R = -6(H 



ttG n ) p 
f 2H 2 ) . 



(37) 
(38) 



If the energy density of the universe is mainly due to 
mat ter (a s in the late accelerated expansion) , we can use 
(16 1 and (14 1 to express p, R, H and H as functions of 



H. Thus we get 



[H 3 gg 12 

i jSHgg' 2 



W 2 g 2 g< + 3H 3 g 2 g"]f"(R) 
R(H) = 18Hg(H)/g'(H) - 12H 2 . 



(39) 
(40) 



H(R) can be obtained from (40 1 and set into (39 1; then 



we get an inhomogeneous second order linear differential 
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equation with non-constant coefficients. Therefore, there 
will always be a bi-parametric family of f(R) actions 
which present ACM as their homogeneous and isotropic 
solution. The difference between these actions will ap- 
pear in the behavior of perturbations, which is not fixed 
by ( 14 ) . Some of these actions are particularly easy to 
solve. If g(H) = H 2 then R = -3H 2 and the differential 
equation becomes 



form of the action will be 



S = 



1 



-g[R + f(G) + WnG N L m ], (45) 



6R 2 f"(R) - Rf'(R) - f(R) + 2R = 0. 



(41) 



from which we derive the modified Friedman equation 

(8irG N ) p=3H 2 - l -Gf'{G) + l -f{G) + l2f"(G)GH 3 . 

(46) 

Following the same procedure as in the case of f(R)- 
theories we arrive to the differential equation 



Its general solution is 

f(R) =R+c x R^ 7 - 



C2 rM?+^) , (42) 



3g(H) = 3H 2 - \Gf{G) + §/(G) 
b^^[9gg' - Hg' 2 - 3Hgg"]f"(G) 

G = 2AH 2 (H 2 - 3Hg/g') . 



(47) 
(48) 



the differential equation (39) will not be solvable ana 



which will give ( 14 1 as First Friedman equation as long as 
the radiation energy densitycan be neglected. In general For a S™E 9(H) , one can use ((48]) to get H(G), then set 

it into ( 47 1 and solve the second order linear differential 
equation. As in the previous case, there will be a bi- 
parametric family of solutions for /(G) which will have 
( 14 1 as their homogeneous isotropic solution as long as 



lytically, and only approximate solutions can be found 
as power series around a certain singular point Rq (a 
value of R such that H(R) cancels out the coefficient 
of f"(R) in ( 39 1 ) . These solutions will be of the form 
f(R) = f P (R) + ci/i(i?) + c 2 f 2 (R) with 



MR) 



oc 



where i 



Ro n 

m— 



p.1,2, 4 1} 



R 



Ro 



-r 



(43) 



matter dominates. The parameters will need to be fixed 
in order to make the contribution of undesired terms in 
the Friedman equation to be negligible when radiation 
dominates. Again g(H) = H 2 is an example which can 
be solved analytically. The differential equation turns to 
be 



,(2) 



1 and the series will 



converge inside a certain radius of convergence. An in- 
teresting choice of R is R Q — —12H 2 _ = R_, which is 
the value of R at H_ . One can also find the approximate 
solution of the differential equation for | R |^>| i?_ |, 
which will be of the form 



12G 2 /"(G)-G/'(G) + /(G) = 0, 



which has the solution 



f(G)= Cl G + c 2 G 1/12 . 



(49) 



(50) 



MR) 

R 



E 

m=0 



^ ){ R- r+m 



(44) 



One can in principle assume that the action (44 1 could 



be considered as valid also in the region in which ra- 
diation begins to dominate, but this action reproduces 
(14 1 only if matter dominates. However, some of the 



In most cases this procedure will not admit an analyti- 
cal solution and an approximate solution will need to be 
found numerically. 

A similar discussion can be made to explain the early 
accelerated expansion (inflation) as a result of an f(R) 
or f(G) action, taking into account that in this case the 
dominant contribution to the energy density is radiation 
instead of matter. The analogue to (44 1 will be now a 
solution of the form 



new terms appearing in (44 1 will be negligible against R 
when radiation begins to dominate. The others can be 
canceled out by setting to zero the appropriate integra- 
tion constant and therefore the Cosmological Standard 
Model will be recovered as a good approximation when 
radiation begins to dominate. 

/(i?)-theories are not the only modified gravity theo- 
ries studied in the literature; /(G)-theories [34] are also a 
popular field of research. In these theories, the Einstein- 
Hilbcrt action is supplemented by a function of the 



MR) 

R 



OC 

n=0 



(— ) 

y R 1 



(51) 



Gauss-Bonnet scalar G = R 2 - AR^R^ 



r> ftV-vpo 
1 *>[iv pa 1 L 



(not to be confused with the gravitational coupling 
Gjv)-The same approach can be used to answer what 
/(G) actions are able to reproduce ACM homogeneous 
evolution, once more in analogy with the modified f(G)- 
gravity reconstruction of a FRW cosmology [35 . The 



for \R\ < \R + \ with R + = -\2H\. Terms in the action 
which dominate over R when R becomes small enough 
should be eliminated in order to recover the Standard 
Cosmological model before matter starts to dominate. 

Both accelerated expansions can be described together 
in the homogeneous limit by an f(R) ac tion with terms 
(^jf) Si+m with Si > coming from (|44| and terms 
(jr-) ri+n with n > coming from (pi|K The action 
will contain a term, the Einstein-Hilbert action R, which 
will be dominant for <C H <C H + including the pe- 
riod when matter and radiation have comparable energy 
densities. 
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ALTERNATIVE DESCRIPTIONS OF THE ACM which, solving for p d , is trivially equivalent to (14 1. This 



In general, by virtue of the gravitational field equa- 
tions, it is always possible to convert a modification in 
the gravitational term of the action to a modification in 
the matter content of the universe. In particular, at the 
homogeneous level, it is possible to convert a generalized 



first Friedman equation of the type ( 14 1 to an equation 



in which, apart from the usual matter term, there is a 
dark energy component with an unusual equation of state 
V — p(p) an d in which General Relativity is not modified. 

The trivial procedure is the following. Assume that the 
source of the gravitational field in the modified gravity 
theory behaves as pd = ujpd (matter or radiation) . This is 
the case when the modification of gravity is relevant. We 
can then define a dark energy or effective gravitational 
energy density as 



Pa 



8irG 



(H 



N 



9(H)) 



(52) 



The continuity equation (11 1 allows to define a pressure 
for this dark energy component as p g — —p g — p g /3H. 
The equation of state of the dominant content of the 
universe leads to 



H 



-3(1 + w) 



HgjH) 
9'{H) 



(53) 



which can be used to express p g as a function of H , and 
one gets the final expression for p g , 



Pg 



-3 



H 2 -2(l+w)H 



9(H) 
9'(H) 



+ ujg(H)\ . (54) 



Then, for any given g(H) we can use (52 1 to get H = 
H(p g ) and substitute in (54 1 to get an expression of 
p g (p g ; H± , to) which can be interpreted as the equation 
of state of a dark energy component. In the simple case 
of a matter dominated universe with ACM, a_ = 1 gives 
obviously a dark energy component verifying p g = — p g . 
Another simple example is a_ = 2, for which 



8ttG 



N 



3H 2 



8ttG 



3H 



^Pg 



(55) 



The inverse procedure is also straightforward. Suppose 
we have a universe filled with a standard component pd = 
ujpd and a dark energy fluid p g — p g {p g ). Using the 
continuity equation of both fluids one can express their 
energy densities as a function of the scale factor a and 
then use this relation to express p g as a function of pd- 
Then the Friedman equation reads 



H 2 



8ttG 



N 



(Pd + Pg(Pd)) , 



(56) 



argument could be applied to reformulate any cosmologi- 
cal model based on a modification of the equation of state 
of the dark energy component |36] as a generalized first 



Friedman equation (14) 



Until now we have considered descriptions in which 
there is an exotic constituent of the universe besides a 
standard component (pressureless matter or radiation). 
In these cases, pressureless matter includes both baryons 
and Dark Matter. However, there are also descriptions 
in which Dark Matter is unified with Dark Energy in a 
single constituent of the universe. One of these examples 
is the Chaplygin gas p g — — l/p g [37j|38]. The use of the 
continuity equation leads to 



Pg 



VA + Ba- & , 



(57) 



where A and B are integration constants. The previous 
method can be used to find the generalized first Friedman 
equation for the baryon density in this model, 



inG N H 2 
— Pb = ~r 




H 4 



(58) 



where H- = y^pA 1 / 4 , k = Bp^ 2 - 1, and p b0 is the 
present value of p\>. In the H ^S> H_ limit this model can 
be interpreted as a universe filled with baryons and dark 
matter or as a universe filled with baryons and with a 
higher effective value of Gat. Equation (58) does not ful- 



fill (22) because it describes the behavior of just baryon 
density. In the H > limit, the model can be inter- 
preted as a universe filled with baryons and a cosmologi- 
cal constant or as an ACM model with a_ = 1 and filled 
only with baryons. 

A similar example in the early period of accelerated ex- 
pansion would be a universe filled with a fluid which be- 
haves as a fluid of ultra-relativistic particles if the energy 
density is low enough but whose density has an upper 
bound 



Px 



C 



(59) 



This dependence for the energy density in the scale factor 
follows from the equation of state 



1 



Px 



;Px 



AC 
3^ 



Px 



(60) 



This model turns out to be equivalent to a universe filled 
with radiation following eq. ( 26 ) with a + = 1 and = 

I S-rGn \ s_ 
\ 3 >C- 

Another equivalent description would be to consider 
that the universe is also filled with some self interact- 
ing scalar field <p which accounts for the discrepancy be- 
tween the standard energy-momentum and GR Einstein 
tensors. Given an arbitrary modified Friedman equation 
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(14 1 a potential V(ip) can be found such that the cos- 



mologies described by both models are the same. The 
procedure is similar to the one used to reconstruct a po- 
tential from a given cosmology From the point of 
view of the scalar field, the cosmology is defined by a set 
of three coupled differential equations 



(p + iH^+V'(ip) = 



IttG 



N 



-3(1 +u)Hp d = p d . 



(61) 

(62) 
(63) 



The solution of these equations for a certain V(<p) will 
give H(t) and Pd(t), and therefore Pd(H) which is g(H) 
up to a factor 87r ^ w . In this way one finds the generalized 
first Friedman equation associated with the introduction 
of a self interacting scalar field. Alternatively, given a 
function g(H) in a generalized first Friedman equation 
(14 1 for the density pd, one can find a scalar field theory 



leading to the same cosmology in the homogeneous limit. 



By considering the time derivative of (62) and using (61 ) 



and ( 63 I one gets 



-(1 + uj)p d 



H 



AttG 



N 



(64) 



where we can use (53) and (14) to get Cp as a function 



of H. Setting this on (62 1 we get V(tp) as a function of 
H. On the other hand, <p'(H) = <p/H, so 



{2H-g'{H))g'{H) 
9(1 + uj)H 2 g(H) 



(65) 



and 
8ttG 



JV 



V(cp(H)) = H 2 -g(H)- 



{l+u){2H-g'{H))g{H) 

2g'{H) 



(66) 

If g(H) is such that the rhs of (65 1 is positive definite, it 



can be solved and the solution ip(H) inverted and substi- 



tuted into (66 1 in order to get V(ip). This is the case of a 
function of the type ( 28 1 with a_ > 1. For the particular 



case of an ACM expansion with q_ = 1, the solution is 
a flat potential V = Vq and a constant value of ip = lpq. 



If g(H) is such that the rhs of (65 1 is negative definite, 
as it happens in ( 28 1 with a_ < 1 or in ( 26 1 , the problem 



can be solved by changing the sign of the kinetic term in 
(62 1. The result is a phantom quintessence in the case of 
(28 1 with a_ < 1. The case (26 1 is more complicated be- 



cause the energy density of the associated inflaton turns 
out to be negative. Moreover, it is of the same order 
as the energy density of ultra-relativistic particles during 
the whole period of accelerated expansion. Therefore, it 
seems that this model is inequivalent to other inflation 
scenarios previously studied. 

In summary, there are many equivalent ways to de- 
scribe the discrepancy between observed matter content 



of the universe and Einstein's General Relativity. At the 
homogeneous level, it is trivial to find relations among 
them. The possibility to establish the equivalence of dif- 
ferent descriptions is not a peculiarity of the description 
of this discrepancy in terms of a generalized first Fried- 
man equation (14 1. The same relations can be found 



among generalized equations of state, scalar-tensor the- 
ories and f(R) modified gravity [40], 



SUMMARY AND DISCUSSION 

It may be interesting to go beyond ACDM in the de- 
scription of the history of the universe in order to identify 
the origin of the two periods of accelerated expansion. 
We have proposed to use the expression of the energy 
density as a function of the Hubble parameter as the 
best candidate to describe the history of the universe. 
In this context the late time period of accelerated ex- 
pansion and the early time inflation period can be easily 
parametrized. 

We have considered a simple modification of the cos- 
mological equations characterized by the appearance of 
an upper and a lower bound on the Hubble expansion 
rate. A better fit of the experimental data can be ob- 
tained with this extended cosmological model as com- 
pared with the ACDM fit. Once more precise data are 
available, it will be possible to identify the behavior of 
the energy density as a function of the Hubble param- 
eter and then look for a theoretical derivation of such 
behavior. 

We plan to continue with a systematic analysis of dif- 
ferent alternatives incorporating the main features of the 
example considered in this work. Although the details of 
the departures from the standard cosmology can change, 
we expect a general pattern of the effects induced by 
the presence of the two bounds on H. We also plan to 
go further, considering the evolution of inhomogeneities 
looking for new consequences of the bounds on H. 

The discussion presented in this work, which is based 
on a new description of the periods of accelerated expan- 
sion of the universe, can open a new way to explore either 
modifications of the theory of gravity or new components 
in the universe homogeneous fluid. Lacking theoretical 
criteria to select among the possible ways to go beyond 
ACDM, we think that the phenomenological approach 
proposed in this work is justified. 
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